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Abstract. In this paper, we present a stability criterion for Processor Sharing 
queues, in which the throughput may depend on the number of customers in 
the system (in such cases such as interferences between the users). Such a 
system is represented by a point measure-valued stochastic recursion keeping 
track of the remaining processing times of the customers. 



1. Introduction 

In this paper, we address the question of stationarity in the general ergodic 
framework for processor sharing queues, in which the throughput (i.e. the quantity 
of work achieved by the server (s) per unit of time) may depend on the state of 
system. More precisely, we assume hereafter that the server(s) (it will be clear in 
the sequel that the effective number of servers does not really matter, only does the 
quantity of work consumed per unit of time) process(es) all the jobs present in the 
system simultaneously and fairly. Whenever there are n customers in the system, 
each of them is thus served at a rate that depend on n, say r(n). The classical 
case is when r(ri) = 1/n, n > 1, so that the total throughput equals n.r(n) = 1 
whenever the system is non-empty: this is the classical Processor Sharing queue. 
Hereafter we consider a more general context, in which the total throughput may 
decrease with the number of customers in the system (hence n.r(n) < 1). This 
is the case for instance in a wireless network in which the number of users being 
currently active may decrease the efficiency of the resources. Another case, is when 
the value of n the number of customers does not change the nominal service rate 
r(n), say r(n) = 1 for all n. This corresponds to the classical queue with infinitely 
many servers. 

In both cases and under general stationary ergodic assumptions, Loynes' stability 
result does not hold since this is not a proper G/G/l queue (the throughput may 
be less, or larger than one). We address the question of the existence of a stationary 
version of such queues by representing them with point measure-valued stochastic 
recursions in the Palm setting, so as to take into account the dependency in the 
number of customers. This point measures keep track of all the remaining service 
times of all the customers in the system. Then, it is possible to provide conditions 
for the existence of a stationary version of this sequence, that allow to explicitly 
construct stationary queues under these assumptions. 

This paper is organized as follows. After some preliminaries in section [21 we 
present the queueing models we consider in section [3j In section [4] we study the 
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particular case of the G/G/oo queue, and in section[5]we present a stability criterion 
for generalized processor queues with state-dependent throughput. 

2. Preliminaries 

Let My and Cb denote respectively the set of positive finite measures on R+ 
and the set of bounded continuous functions from K to K.. Equipped with the weak 
topology a ^My ,Cbj , My is Polish (see [2]). Let be the zero measure on K (i.e., 

such that 0(23) = for any Borel set 55 on R). For any [i S My and any measurable 
/ : R — > R, we classically write (/i, /) :— J f d[i. Let us denote for any y G M and 
any measurable / : R — > R, T y f(.) = /(. — y)l{ >y }. Then, for any /i S Mi, 
denotes the only element of My s.t. (r a ^, /) = (^,T y f). 

Let the set My be endowed with the increasing partial integral order ^ : for any 
two [i, v G My , ^ if /) < (i/, /) for any measurable non-decreasing function 
/ such that these integrals exist. Of course, < /i for any ^ S My . Furthermore, 
let us remark that 

Lemma 1. Any sequence o/My that is ^-increasing and bounded above converges 
for the weak topology. 

Proof. Let {/^nlngpj be a ^-increasing sequence of My that is bounded above 
by /i £ My . Then, as easily seen the sequence of non-increasing real functions 
{\i n ([-,oo))} neN tends pointwise, and hence (this is Diniz Theorem), uniformly to 
a non-increasing real function / that is right continuous and has a countable num- 
ber of discontinuities. Moreover /(0) < /i(R!j_) < oo, and we can fully characterize 
a measure fi* G My setting fi* ((0, x)) = /(0) - f(x) for all x G R+. In particu- 
lar, sup^gR. \fj, n ((0, x)) — /i* ((0, a;)) I — > 0, hence /i™ tends to fi* in total variation, 

+ n — >oo 

which completes the proof. □ 
Let now M. C My be the subset of finite (simple) counting measures on W + . Any 

/i G A / (\{0} reads [i = J2i=i ^0)' wnere :— /i(R+) is the number of atoms 

of /i, S x is the Dirac measure at x G R+ and ai(/i) < 02 (aO < ■•• < o^N^if 1 )- Then, 
T vip) — SfcLi 3ai{ti)-y~L{ai(p)>y} and for any two fi, v G Ai\ {0}, [i < v whenever 

(i) N(n) < N(y), 

(ii) for all i = 0, ...,N(n) - 1, ajv( M )-i()") < 

We denote for any G {0}, Z{\i) = a A r( At )(/x), the largest atom of /x. Finally, 
we write x + = max(x, 0) for any real number x, 'Yin— a ■ = whenever k < j and 
max{0} = 0. 

3. The model 

Let us first introduce our definitions and assumptions on the queueing systems 
we shall consider in the sequel. Let (O, P, 9 t ) be a probability space furnished 
with a bijective flow (9t) t>0 , under which P is stationary and ergodic. Define on Q 
the ^-compatible simple point process (A t ) teR of points ... < T_2 < T-\ < Tq < 
< T\ < T2 < that represent the arrival times of the customers in a queue 
without buffer. The process (A t ) t£R is marked by a sequence {cr n } n ^ where for 
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all n G Z, a n is the service duration requested by the customer C n arrived at time 
T n . Also denote for all n G Z, £„ = T n+ i — T„, and suppose that the generic r.v. a 
and £ are integrable. We consider that the server(s) follow a generalized Processor 
Sharing discipline. By that, we mean that all present customers are taken care of 
simultaneously, at a rate r that is equal for all customers. An example is of course 
provided by the classical Processor Sharing queue, but it will be shown in the 
subsequent sections that significant results can be obtained as well for a wider class 
of systems. Indeed, it is plausible to assume in many cases, that the amount of work 
in the system might affect the throughput, considering for instance the working cost 
induced by the switching mechanism in the processor, or the interferences between 
the users of a wireless network. In both cases, it is then natural to assume that 
the rate r is a non- increasing function of the service profile, i.e. [i < v implies 
r (p) > r i v )- Hereafter, for the sake of simplicity, we will restrict to the sub-case, 
where r is a non-increasing function of the number of customers in the system, 
although it should be clear that all the results below hold as well when r is function 
of the whole service profile. In other words, at any t, denoting Q{t) the number of 
customers in the system at t, each customer is allocated a quantity of work r(Q t ) 
per unit of time, that is such that r(i) > r(j) for all i, j G N* such that i < j. Let 
us illustrate through a naive example the effect of a large number of customers on 
the throughput. 

Nominal service rate Throughput 

1 customer 1 1 

2 customer 0.495 0.99 

3 customer 0.3 0.9 

100 customers 0.008 0.8 

Provided that C n is in the system at t, his remaining processing time at this 
instant is the time before his service completion. The service profile of the system 
at t is the A4-valued process keeping track of the remaining processing times of all 
the customers in the system at t: 

Q(t) 

M*) = ^2 <*«i(M(t)) 

i=l 

where a\{ii{t)) < a2(^(t)) < < aQ^(^i(t)) denote the remaining processing 

times of the customers in the system at t, ranked in decreasing order. Let W(t) 
denote the workload at t. Then, the workload and the congestion processes can 
easily be recovered from the service profile process by writing for all t 

Q(t) =N(n(t)), 
W(t) =(n{t),I), 

where I is the identity function. The processes n, Q and W have RCLL paths, and 
we denote for all t n{t—) = lim s || t /j,(s) (and accordingly, Q(t—) and W(t— )). We 
denote for all n G N, fj, n = /x(T n -) (respectively Q n = Q(T n -), W n = X(T n -)) 
the service profile (resp. congestion, workload) just before the arrival of customer 

Let (yijJ 7 , P°) be the Palm space of A, denote 8 := 8t x , 0" 1 his measurable 
inverse and for all n G Z, 6 n = 6 o 6 o ... o 6 and 6~ n = 0- 1 o Q- 1 o ... o Note, 
that P° is stationary and ergodic under 9, i.e. for all 21 G P° [9- 1 ^i] = P° [21] 
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and 6*21 = 21 implies P° [21] — or 1, and that all ^-contracting event (such that 
P° [2l c n 6'~ 1 2l] = 0) is ^-invariant. Denoting £ := £o and a := Co, we have for all 
n G Z, £ n := £ ° and a n := a o 0". 

We say that the £7- valued random sequence {X n } ngN is a stochastically recursive 
sequence (SRS) whenever for some random mapping <fi '■ E — > E, 

X n+1 =<f>o6 n (X n ) , n G N, P° - a.s.. 

For any £7-valued r.v. F, we denote { X„ } the SRS {X„} n eN such that X^ = 

y , P -a.s.. We follow the formalism of [1] and formulate the question of stationarity 
for the SRS {^ n } TlgN in the following terms. There exists a stationary version of 

{X n } n&N whenever for some Y and for all n, X l J ] = Y o n , P° -a.s., or in other 
words, provided that the equation 

Y o0 = (f)(Y) 

admits a solution that is a E- valued r.v.. We say that two sequences of r.v. {X n } neN 
and {Y n } n€N couple provided that 

P° [BN(w), X n {u) = Y n {u) for all n > N(lu)} = 1, 

and that there is strong backwards from {X n } nefi with the stationary sequence 
{Y o 6 n } whenever 

P° [3N'(w), X n o e- n (u) = Y{uj) for all n > N'(lu)] = 1. 

Lemma 2. The sequence {^ n } ne ^ is stochastically recursive for any rate function 
r: letting for all /i € M. and x G R+, 
• For all i < N(n), 

Yii^x) = r(JV(/i) - * + 1) fx- ^a^Oi) 
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K r(N(ri-j + l) r(N(ji)-j) 

• i r (fi,x) = maxji < ATQu); a^n) < 7[(/x,a;) j, 

• 7 r (A^) ==7(iro,, x )+i)AiO*,aO, 

we have for any initial profile /j,q anf for all n G N, 

(1) Mn+l = * r (Mn + <W , £n) • 

Proof. Just after the arrival of C„, the service profile reads \x := fx n + 8 rTn . Set 
Tq := T n and a (ju) = 0. For any % G {1, iV (/*)}, let T[ be the theoretical 
departure of the customer (7, whose remaining service time at T n is a,(/i). The 
remaining service time of at T[_ Y is a,;(/i) — aj_i(/i), and between T[_ Y and T/, 
Ci is served at rate r(N({i) + Hence we have the induction formula 

(2) Ti=n- 1+ a f^- ai ;^ ,i^{i,..,Nm, 

r(N(fi) - i + 1) 
from which we deduce that for all i € {1, N(fi)}, 



( 3 ) n=T n+ ";M^ + Y,«M 



i-l 



r(N(fi)-i + l) ,w W%)-jtl) r(N(jj)-j) / 
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For any i, customer G% leaves the system before T n+ i provided that T[ — T n < £„, 
which is equivalent to Oi(fi) < 7[(/i, £ n ) m view of ([3]). In particular, « r (£t,£n) 
denotes the index of the last customer leaving the system before T„+i (or if there 
is no departure between T n and T n+ i). Then the system is not empty at T n+ i — 

provided that z r (^,£„) < N(fj,), and in this case, (Ci,i € {i r (^,£n) + Ij-^G 14 )}} is 
the set of customers present in the system at T n+ i — . For such i > i r (fi,£ n ), the 
remaining service time of Ci at T n+ i is given by 

oti(n) - atr (AtiCn )(^) - r(N(n) - i r (/i,£„)) (r n+1 - ^r (/lifn) J = a^fi) - Y (//,£„) . 
Thus the functional mapping the profile at T n onto the profile at T n+ \— reads 

N(p) 

$ r (.,£„) : /i i — * ^2 <Wp)-7 r (Mn)- 
i=i r (M,6>)+l 

To obtain the announced result, remark that for any S M. and a; S R+, for any 
i < AT(/x) we have that 

r(N(fi)-i)-r(N(v)-i + l) , r{ \ ( ,\\ 
Jt+iiv-.x) -7i(M.«) = rTM/i) - i + 1) ('W' 35 ) _Q! i(M)), 

which is nonnegative if and only if i < i' r (/J, x). Hence, 
(4) Y (ji, x) = max 7 [ (/x, x) , 

l<i<Af(jii) 

and in particular $ r (/!,£„) = ^(^^jM, P°-a.s.. □ 

For a fixed x e M+, the two following monotonicity properties of the mappings 
<I> r (.,x) hold, as shown in Appendix. 

Lemma 3. For any x G R+ and any rate function r, the mapping <I> r (.,x) is 
non- decreasing from Ai into itself. 

Lemma 4. For any x G R+ and any /j, £ M, for any two rate functions r and f 
such that r(i) < f(i) for all i £ N*, ^ r (^, x) h 3> f (/J, x). 



4. The pure delay system 

Let us first consider the case, where the rate function is constant with respect to 
the size of the system, say r(i) — 1 for any i > 1. This corresponds to the classical 
"pure delay" G/G/oo queue: all present customers are simultaneously served at 
unit rate, and hence spend in the system a time equal to their service duration, 
which is equivalent to say that there is an infinity a servers. In this case, the 
recursive equation <[TJ) driving the service profile sequence (for which a diffusion 
approximation is given in [5] in the M/GI/oo case) specializes to 

(5) Mn+l = T «„ (j"n + 5<r n ) 

and a stationary service profile for the queue is a solution to the equation 

(6) fl o 9 — 7-e (fJt + 8a) . 

The following lemma (see [7]) will be used in the sequel. 
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Lemma 5. There exists a unique P°-a.s. finite solution to the equation 
(7) Lo9= [max {L, a}- £] + , 

given by 



(8) 



L := 



sup a_j - £-i 



Proof. Existence. Loynes' Theorem for stochastic recurrences (see [6], [T]) can be 
applied since the mapping x i— > [max{x, ct} — £] + is P°-a.s. continuous and non- 
decreasing. The minimal solution L to ([7| classically reads as the P°-a.s. limit of 

Loynes's sequence 1 Lj? r , where < > is the initially null SRS that 

is defined by 



L 



[o] _ 

n+1 — 



for all n £ N. 



It is routine to check from Birkhoff 's ergodic theorem (and the fact that a is not 
identically zero) that L is P°-a.s. finite. 

Uniqueness. Let L be a solution to (7j). First, remark that if L > a, P°-a.s. would 
imply that on a P°-a.s. event, 



Io0>OoIo( 



a contradiction to the ergodic Lemma. Hence in view of the minimality of L, we 
have that 



L = L 



> P L 



L < a 



>0, 



P u Lo6 <Lo, 

which implies that |i = L,\ is P°-almost sure since it is ^-contracting. 

We can now state the following result. 
Theorem 1. The equation admits a finite solution, given by 



□ 



i=l 

Moreover, provided that 

(9) P° [L < 0] > 0, 

this solution is unique and for all £ such that Z(() < L, P°-a.s, the sequence 
luj? \ converges with strong backwards coupling to [i PD . 

Proof. Existence. It is a straightforward consequence of Birkhoff 's ergodic theorem 
that 



P° [/x PD 6 M] 



Card ^ t e N*, a-i - £-i > f < 

3=1 



> 0. 



This 0-contracting event is thus P°-almost sure. On another hand, in view of 
Lemma O the mapping fx i— > rg(/i + ^ CT ) is P°-a.s. non-decreasing from into 
itself. It is furthermore continuous for the weak topology, as easily checked from 
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the fact that for any ./W-valued sequence {^n}„ eN tending weakly to v, for any 
x, s e R+ and any cj) £ Cb, 

{t x v„ + S s ,<(>) = I 4>{y - x) dv n {y) + <p(s) — > / <f>(y - x) dv(y) = (r x v + 5 s ,(p). 
J n ^°° J 

Thus, we can follow the steps of Loynes' construction (Lemma [T]), to conclude that 
/x PD is the ^-minimal solution of <(Gj since it is the P°-a.s. limit of the sequence 
given for all n e N by 

oo 

Uniqueness. It is easily checked that for any solution /i of ©, 

Z(ji) o9 = Z( Ti (n + S a )) = [Z{ji) V a - £]+ , 

hence Z{p) = L, P°-a.s.. Moreover, since ^ PD is the minimal solution of (JH), we 
have that 

V = M PD } _> {M = 0} - {Z{n) = 0} = {L = 0} . 
Hence, whenever <(9j) holds, the event {/i = ^t PD } has a positive probability. Since 
it is ^-invariant, it is P°-almost sure. 

Coupling. Let £ be a .M-valued r.v. such that Z(Q < L, P°-a.s.. It is easy to 
construct another .M-valued r.v. C such that £ C an d Z(() = L, P°-a.s. by 
setting e.g. £ = J^fJP 1 ^»(C) + From Lemma it follows by induction that 
mL CI d 1$, P°-a.s. for all neN. Remark now that for all n 6 N, Z (/4£ ] ) =Lo9 n , 
as easily checked by induction. Hence, for all n G N, we have 

S n := {L o <r = 0} = {z (^p) = o} = {/4? = o} C {/4« = o} . 

Therefore, {£n} ne pi is a stationary sequence of renovating events of length 1 for 
(see Oil]) for any C such that Z(Q < L, P°-a.s.. Assumptions IM 

L J neN 

implies the coupling property for such an initial condition in view of Corollary 
2.5.1 of pQ. □ 

As simple consequences of the latter result, let us remark the following coupling 
properties. 

Corollary 1. Under condition (G|), for any £ such that Z{Q < L, P°-a.s, 
(i) <Jfn } converges with strong backwards coupling to N {n PD ) ; 

I J neN 

(a) \wk {c ' I)] \ converges with strong backwards coupling to (/i PD ,I). 

I J neN 

5. Processor Sharing queues 

We shall now consider the case, where the rate function depends on the number 
of customers in the system at current time. We assume hereafter that the non- 
decreasing function r is such that 

(10) sup n.r(n) < 1, 

neN* 

(11) K r = inf n.r(n) > 0. 

neN* 
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Assumption (flO|) amounts to say that there is a single server, in that the throughput 
at time t, given by Q(t).r(Q(t)), may not exceed one. A typical case is the classical 
Processor Sharing queue: assume that r(n) = rT x for any n (and hence K r = 1), 
meaning that all customers are served at a rate that is inversely proportional to 
the number of customers. In that case the server works at unit rate whatever the 
number of customers in the system. Whenever K r < 1, the number of customers 
affects the velocity of service, so that the total throughput may be less than one. 
We assume nevertheless in l[TTj) that a minimal throughput K r is granted for a 
given r, i.e. the server always achieves at least K r unit of work per unit of time. 
An example is provided by the following idealistic scenario: the server works at 
unit rate whenever there is only one customer in the system (r(l) = 1), and the 
interferences (or operating cost) when there are several customers in service at the 
same time decreases by half the efficiency of the server, so that r(i) — 1 / (2i) for 
any i > 2, which implies in particular that ifTTj) is met for K r = 1/2. 

In view of Lemma [2j a stationary service profile is a solution to the equation 

(12) tioe = $ r (n + 6 a ,0- 

We have the following result. 

Theorem 2. Let r be a rate function satisfying assumptions flO\ ) and (Jij) . Then 
provided that 

(13) E° [a] < K r E° [£] , 

the equation \12jl admits a unique finite solution fx r . Moreover, for any Ai-valued 
r.v. ( such that ((,1) < W r , P°-a.s. (where W r is the unique solution of 

the sequence \un \ converges with strong backward coupling to if . 

I J riGN 

Proof. Existence. Fix r satisfying (fTUj) and (TTJ . From Loynes's fundamental sta- 
bility result, the equation 

(14) Wo6=[W + a- K r S] + 

admits a unique P°-a.s. finite solution, say W Kr , provided that lfT3|) holds. Let f be 
the rate function such that for all fi € M, f(p) = K r /N{fj), so that the throughput 
under f always equals K r whenever the system is non-empty. Let C be a ./Vf-valued 
r.v. such that ((, I) < W Kr and 

C = C + ^w k ^~{c,i)^-w k ^>{cj)- 
Is then clear that (C, I) = W Kr . Moreover, we have P°-a.s. for all n £ N 

as the throughput equals K r at any time (as easily checked from Lemma [2|, so that 

(Mn+i>-0 = W Kr o 9 n for all n £ N. On another hand, C < C, hence in view of 

Lemmas [3] and [H an immediate induction shows that Mn [C1 d M« [C1 for all n e N, 
which implies in turn that 

(A# [<] > J) < (»n®,I) = W Kr o6 n for all n € N. 
Therefore, for all n G N, on 2l„ := (W^ o 6>" = 0}, we have that (^' [C1 , 1) = 0, 

r \0 

hence // n = and 
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Therefore {/z„' lCI } admits {2l„}„ eN as a stationary sequence of renovating events 

L J neN 

of length 1. Furthermore, the event 2lo = {W Kr = 0} has a striclty positive prob- 
ability since the contrary would imply that 

E° [W Kr 06- W Kr ] = E° [a - K r £] < 0, 

an absurdity in view of the ergodic Lemma. Then it follows from pQ, Th. 2.5.3, 

r [CI 

that there is strong backwards coupling of jU„ n with the stationary sequence 



{/x r o # n } ngN , where [f is a proper solution to lfT2|) . 

Uniqueness. Fix r and f be as above. There exists a solution [i r to (fl2"ll . Then, 
we have P°-a.s. 

(//, 7) o = <<F (// + S a , £) , I) = [<//, /) + a - KrZ] + , 

hence I) equals W Kr , P°-a.s.. Moreover, on {{fi r , I) < VF^ } , we have in view 
of Lemma [2] that 

(ff,I)°0< (<F (// + 8„, , 1) = [iff, I) +a- K r (\ + < W K " a 6, P° - a.s., 

thus the event {(/i r ,i) < W Kr } is ^-contracting. Moreover, 

P° [<//,/) <W^] >P° [(^,I) = 0}>0, 

as another consequence of (fl3"l) and the ergodic Lemma. Therefore, (fJ. r ,I) < W Kr , 
P°-a.s., so that 

On C {< M r , /> o 0" = 0} = {// o 0" = 0} . 

Consequently, {2ln}„ eN is a stationary sequence of renovating events of length 1 
for {/i r o #™} neN for any solution fi r of the equation lfT2|) associated to the rate r. 
Since P° [2lo], there exists a unique solution to lfT2|) in view of Remark 2.5.3. in 

ED- □ 

We have in particular: 
Corollary 2. Under condition for any C such that (£,1) < W Kr , P°-a.s., 
converges with strong backwards coupling to N (/i r ) ; 

I J neN 

(a) {wi {u)] \ converges with strong backwards coupling to 

Appendix A. Proofs of monotonicity 

For easy checking, we present hereafter the details of the derivations proving 
Lemmas [3] and [H 

Proof of Lemma\^ We fix again x 6 R+ and fi, v E M such that [i<v. Whenever 
i r {^ji,x) < N(fi) (otherwise $ r (/i,x) = 0), we have that 
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which implies that 

>- r<NW - (' - 'I? "'<"> U(,)-,- + i) - ?W5P7))) 

^ 77V(i/)-Af( Al )+i''( Al ^) + l( Jy ' 

This means that i (v, x) < N(v)-N((i)+i (n,x), i.e. N ($ r (jU, x)) < N($ r (v,x)). 
Hence in view of ([4]) , we have 

70*i = 7^,0+1 (P,x) > Y i:l r {L ,, i)+N(p) „ NH)++1 {v,x) 

= 7 r (^, b), 

which clearly implies that < & r (^,a;) X $ r (V, x). □ 

Proof of Lemma^ We now fajieM and a; £ K+. For any two rate functions r 
and r such that r(i) < r(i) for any ieF, the induction formula ([21) straightfor- 
wardly shows that i r (pL,x) > i r {^,x) i.e. N ( < & r (/i, x)) < N (<I> r (/x, x)) . Hence, as 
in the previous proof, 7 r (A*, x) < 7 r (/U, x). □ 
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